home *** CD-ROM | disk | FTP | other *** search
/ WINMX Assorted Textfiles / Ebooks.tar / Text - Mathematics - Numerical Mathematics and Computing (F).zip / string.f < prev    next >
Text File  |  2002-06-11  |  1KB  |  47 lines
  1. C
  2. C PAGE 441-442: NUMERICAL MATHEMATICS AND COMPUTING, CHENEY/KINCAID, 1985
  3. C
  4. C FILE: STRING.FOR
  5. C
  6. C HYPERBOLIC PDE PROBLEM SOLVED BY DISCRETIZATION (F,TRUE)
  7. C
  8.       DIMENSION  U(11),V(11),W(11)    
  9.       DATA  N/11/, M/20/, H/0.1/, HK/0.05/      
  10.       DATA  U(1),V(1),W(1),U(11),V(11),W(11)/6*0.0/       
  11.       RHO = (HK/H)**2       
  12.       PHO = 2.0*(1.0 - RHO) 
  13.       DO 2 I = 2,N-1
  14.         X = REAL(I-1)*H     
  15.         W(I) = F(X) 
  16.         V(I) = 0.5*( RHO*(F(X-H) + F(X+H)) + PHO*F(X) )   
  17.    2  CONTINUE
  18.       DO 5 K = 2,M
  19.         DO 3 I = 2,N-1      
  20.           U(I) = RHO*(V(I+1) + V(I-1)) + PHO*V(I) - W(I)  
  21.    3    CONTINUE
  22.         PRINT 6,K,(U(I),I = 1,N)      
  23.         DO 4 I = 2,N-1      
  24.           W(I) = V(I)       
  25.           V(I) = U(I)       
  26.           T = REAL(K)*HK    
  27.           X = REAL(I-1)*H   
  28.           U(I) = TRUE(X,T) - V(I)     
  29.    4    CONTINUE
  30.         PRINT 6,K,(U(I),I = 1,N)      
  31.    5  CONTINUE
  32.    6  FORMAT(//5X,I5,//(4(5X,E22.14)))
  33.       STOP
  34.       END 
  35.         
  36.       FUNCTION F(X) 
  37.       DATA  PI/3.14159 26535 898/     
  38.       F = SIN(PI*X) 
  39.       RETURN      
  40.       END 
  41.         
  42.       FUNCTION TRUE(X,T)    
  43.       DATA  PI/3.14159 26535 898/     
  44.       TRUE = SIN(PI*X)*COS(PI*T)      
  45.       RETURN      
  46.       END 
  47.